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ABSTRACT 


An improved model is developed for predicting the shape of an electric arc balanced 
in transverse convective and magnetic fields. The improvement over the previous model 
is in the more accurate description of the external flow field. The techniques used are 
based on conformal mapping and are an extension of those methods developed for prob- 
lems of free surface flows. 


SHAPE OF A MAGNETICALLY BALANCED ARC 
by Marvin E. Goldstein and Willis H. Braun 
Lewis Research Center 

SUMMARY 

An improved model is developed for predicting the shape of an electric arc balanced 
in transverse convective and magnetic fields. The improvement over the previous model 
is in the more accurate description of the external flow field. The techniques used are 
based on conformal mapping and are an extension of those methods developed for prob- 
lems of free surface flows. 


INTRODUCTION 

The importance of electric arcs interacting with transverse convective and magnetic 
fields has greatly increased during the last decade because of numerous applications of 
electric arc devices to aerospace technology (ref. 1). The two types of configurations 
that have been used to study this threefold interaction are (1) the so-called rail accelera- 
tor in which the arc moves through a quiescent fluid and (2) the magnetically balanced arc 
that is discussed herein (see fig. 1). In this type, a magnetic field perpendicular to both 
a channel flow and an electric field is adjusted so that the electromagnetic force j x B 
acting on the arc just balances the hydrodynamic forces. The experiments show that this 
type of arc tends to behave much like a cylindrical solid body (ref. 1). 

A simple model, from which the cross-sectional shape of the arc could be deter- 
mined, was proposed in reference 2. That paper pointed out that some of the features of 
the model were expected to be inaccurate in detail but that the type of analysis employed 
was expected to be correct. In the present report, the arc model described in refer- 
ence 2 is improved by eliminating some of the inaccuracies concerned with the external 
flow about the arc, although the same model for the internal structure of the arc is re- 
tained. 

Those features of the model presented in reference 2 that are used herein are the 
assumptions that the hot fluid recirculates within the arc column while the cold external 
fluid flows around the column, much as it would flow about a solid body; and, further- 




more, the arc can be considered two-dimensional (all variations occur in the x, y-plane 
shown in fig. 1); the current density j is uniform across the arc cross section, and Hall 
currents can be neglected. Under these conditions, the magnetic force per unit volume 
is constant and equal to ijB, where i is the unit vector in the direction opposite to the 

a 

upstream velocity (as shown in fig. 1). Now ijB is a conservative force and, hence, can 
be balanced by a pressure gradient. Reference 2 also assumed that both inertial and 
viscous forces could be neglected within the arc. With these approximations, the mo- 
mentum equation for the hot fluid within the arc is 


Vp. = ijB (1) 

where p^ is the pressure within the arc. The external flow was taken to be inviscid and 
incompressible. (Most experiments were carried out at sufficiently low velocities so that 
compressibility effects were not important). The purpose of the present analysis is to 
find the shape of the arc that will make the external potential flow compatible with the so- 
lutions to equation (1) within the arc. (The matching conditions are discussed in the next 
section.) At this point, it is desirable to deviate from the analysis given in reference 2. 
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The model for the potential flow given in reference 2, although it satisfied the appropri- 
ate matching conditions along the surface of the arc, had several serious drawbacks. 
First, the arc turned out to be a semi- infinite body that had to be artificially truncated. 
Second, the streamlines gradually converged at large distances from the arc instead of 
becoming parallel. Third, the model did not take into account the wake region behind the 
arc. Fourth, the model differed from the experimental situation in that it did not take 
into account the presence of channel walls that could have a large effect on the wake re- 
gion, even if the walls are quite far apart. In the present model, all these drawbacks are 
eliminated. 


SYMBOLS 


A cross-sectional area of arc 

B |B| 

B magnetic flux density vector 

D channel width 

d longest dimension of arc in direction parallel to flow 

h distance between stagnation point and point at upstream infinity in t- plane 

J total current in arc 

i unit vector in direction opposite to upstream velocity 

i HI 

J arc current density vector 

l half the longest dimension of arc in direction perpendicular to flow 

p pressure in external flow 

p. pressure in arc 

p Q stagnation pressure 

s Re (t) 

s* distance from stagnation point measured along boundary streamline of arc 
T intermediate complex variable 

Tj complex coordinate of point I in T-plane 
t intermediate complex variable 

u x-component of velocity of external flow 
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magnitude of velocity of external flow, | 

magnitude of velocity of external flow on boundary streamline of wake 
magnitude of velocity of external flow on boundary streamline of arc, |?(s) | 
magnitude of velocity far upstream of arc, |?(ih) | 
y-component of velocity of external flow 
complex potential, $ + i^ 

coordinate perpendicular to flow and parallel to magnetic field 

x-coordinate of boundary streamline of wake 

x-coordinate of boundary streamline of arc 

coordinate in direction opposite to flow 

y-coordinate of boundary streamline of wake 

y-coordinate of boundary streamline of arc 

complex variable, x + iy 

x f + iy f 

x s + *8 

included angle between two tangents to arc surface at stagnation point divided by n 
interior of first quadrant of t- plane 
convergence factor 

complex conjugate velocity, u - iv = dW /dz 


1 + 1 T t - 1 


t 

Re (SI) 

positive real number 
density of external flow 
Im (SI) 

velocity potential 

angle between tangent to surface of arc and negative y-axis 
stream function 
V D/27rh 2 
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fi i In (?/V 0 77) 

~ asymptotically proportional to 

Superscripts: 

— arc 

* complex conjugate 


ANALYSIS 

Boundary Conditions on Velocity 

Under the assumptions just discussed and from the integration of equation (1), the 
pressure p^ at any point with coordinates x, y (fig. 1) lying within the arc is given by 

Pi - P Q = iBy (2) 

that is, the pressure in the arc is independent of x. This varying pressure in the arc 
must be balanced at the surface of the arc by the pressure field created by the external 
flow, and the arc surface must correspond to a streamline of the external flow. The dis- 
tance measured along this bounding streamline (the arcs BA or BC in fig. 2) is denoted 
by s* and the magnitude of the velocity of the external flow along this streamline by V g . 
Since the external flow was assumed to be incompressible and inviscid, Euler's equation 
along this streamline becomes 


-pV s 


9V 


s _ 9p 


9s 1 


9s’ 


(3) 


Now, since the internal and external pressure must balance along this bounding 
streamline, 


P = PjCy) on BC 


(4) 
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Figure 2. - Z-plane. 


so that using equation (2) results in 

_§ p _ = = *l 

3s 1 " 3s 1 " d y 3s 1 " 

= jB = -jB cos <p (5) 

3s 1 " 

where <p is the angle between the tangent to the boundary streamline and the negative 
y-axis as shown in figure 2. Substituting equation (5) into equation (3) yields 

pV_ dV = jB cos <p ds 1 ^ (6) 

s s 

where dV is understood to be the change in V 0 along s 1 ". 
s s 

Behind the arc is a wake. The temperature in the arc is very high. The fluid tem- 
perature in the wake is also high, but is somewhat cooler than that in the arc. Since the 
electrical conductivity of the fluids of interest for arc experiments (argon, air, etc. ) 
drops rapidly to zero at a certain temperature and since the temperature must eventually 
drop off behind the arc, there will be a line on one side of which the fluid will be conduct- 
ing and on the other side of which it will be nonconducting. This line is defined as the 
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front boundary of the wake. Its location will be determined by the calculation; it is only 
necessary to assume that such a line exists. Thus the fluid in the wake is nonconducting 
so that the pressure acting on it is caused by only hydrodynamic forces. But this fluid is 
very hot (ref. 1), and therefore its density and the inertial forces are very low. Hence, 
to a very good approximation, the pressure is constant in the wake. The bounding stream- 
lines of this wake must begin at the end points of the arc boundary C and A (i. e. , they 
must be continuations of the bounding streamlines of the arc) and extend to -°°, as shown 
in figure 2. Along these streamlines, the pressure on both sides must balance so that the 
pressure in the external flow is a constant along the streamlines CH and AH bounding 
the wake. Therefore, from Euler’s equation, the magnitude of the velocity of the external 
flow is constant along these streamlines. This constant velocity is denoted by V Q , and 
the arc is assumed to be situated in a channel whose walls are separated by the distance 
D. In the final formulas, taking the limit D — °° will include the case where the arc is 
situated in an unbounded parallel flow. The boundary condition for the velocity along the 
channel walls is determined from the fact that it must be in the -y-direction, as must the 
velocity along the stagnation streamline IB (c. f. , fig. 2). 


Boundary Conditions on Stream Function 

In order to completely determine the problem, it is sufficient to specify the values 
of the stream function it along the boundaries of the flow or equivalently along the bound- 
aries of half the flow, since it is symmetric with respect to the y-axis. To specify these 
values requires that the bounding streamline of the arc BC be a continuation of the stag- 
nation streamline IB and that the bounding streamline of the wake CH be a continuation 
of the arc BC. Thus, ^ = 0 on IB, BC, and CH. In addition, the channel wall corre- 
sponds to a streamline, or St is constant along the channel wall IH. 


Calculation of Stream Function 

The problem is now completely specified, and its solution is obtained by conformally 
mapping the complex potential W defined by 

W = t + rt (7) 

where <t is the velocity potential, and the complex conjugate velocity ? which is defined by 

?= dW = u-iv (8) 

dz 

where z = x + iy, and u and v are the x- and y- components of the velocity, respec- 
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tively, into the upper half of an intermediate t-plane to be specified next. 

Once W and £ are known as a function of the complex variable t, z can be com- 
puted as a function of t from 


z = 



dz dW _ 
dW dt 



r^dt 

dt 


( 9 ) 


Thus, at least in principle, if t is eliminated, W and £ can be found as functions of z. 
This yields the solution of the problem. 

First consider the W- plane. From the previous discussion, the arcs IB, BC, and 

/-'-v 

CH must lie along the line S? = 0, and the wall streamline IH must lie along a line 
’Jr = constant. Thus, as shown in figure 3, the interior of half the flow field must be map- 
ped into an infinite strip in the W-plane. Now, the t-plane is chosen so that the bound- 
aries of the flow field are as shown in figure 4 ; that is, die bounding surface of the arc is 





‘ lm(t> 
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Figure 4. 

- 1- Plane. 
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mapped into the interval (-1, 1) on the real axis; the bounding streamlines of the wake 
are mapped into the rest of the real axis, the channel walls are mapped into the portion 
of the imaginary axis above the point I; and the stagnation streamline is mapped into the 
rest of the positive imaginary axis. The point at upstream infinity in the physical z-plane 
is mapped into the point I on the imaginary axis, and the point far downstream in the 
physical plane is mapped into the point at infinity in the t-plane. 

A mapping must be found of the interior of the infinite strip in the W-plane into the 
first quadrant of the t-plane. To accomplish this maping, note that the transformation 

T = t 2 (10) 

maps the first quadrant of the t-plane into the upper-half T-plane, as shown in figure 5, 
so that the positive imaginary axis in the t-plane is mapped into the negative real axis in 
the T-plane. Now, map the interior of the infinite strip in the W-plane into the upper- 
half T-plane by the degenerate Schwarz -Christoffel transformation (ref. 3, p. 226) 



Eliminating T between equations (10) and (11) and noting from figure 4 that, at the point 
I in the t-plane, t = ih yield 


dW _ 

2 V 


dt 

t 2 

1 + — 
h 2 



i 1m (T) 


I 

B C 

Re (T) 


h 2 —| 0 


( 12 ) 


Figure 5. - T-plane. 
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To establish the connection between h, and the physical variables, integrate equa- 
tion ( 12 ) to get 


W = h 2 * 0 In (t 2 + h 2 ) + constant (13) 

Now, since W(-t*) = W*(t), it is clear from the symmetry of the problem that equa- 
tion (13) gives the correct complex potential for the entire flow field when it is extended 

from the first quadrant of the t-plane to the upper-half t-plane. The imaginary part of 

2 

W can be shown to change by 27rh in passing from the negative to positive side of the 
imaginary axis with Im (t) > h; that is, the change in ^ between the channel walls is 
27rh 2 But this is the total flow through the channel and is therefore equal to V^D, 
where is the upstream velocity of the flow. Hence, it follows that 

VJ) = 2rrh\ (14) 

Equation (12) shows that, in the limit h — °°, 


w = v 2 

which maps the entire W-plane into the upper-half t-plane. This expression is the proper 
transformation of the W-plane into the t-plane for an unbounded flow (channel walls infi- 
nitely far apart). Note that in this case the stagnation streamline maps into the entire 
positive imaginary axis. Hence, at least as far as W is concerned, the limit D — °° 
corresponds to the limit h — The next section makes apparent the fact that £ de- 

pends on h only through dW/dt and so the same conclusion follows for £. 


Calculation of Velocity 


With W 
function of t. 


specified as a function of the intermediate variable t, £ is now found as a 
This is accomplished by studying the properties of the function 77 , defined by 


77 (t) = 



e i(l- 0 )jr /2 


0 1 


(15) 
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The branch of yt^ - 1 is chosen so that the cut lies along the real axis between - 1 
(see ref. 4, p. 82). When t approaches the real axis through values with positive 
imaginary part, 


2 1 1 / 2 

i 1-r -1 < t < 1 


it 2 -ii 1/2 t>x 


i 2 i l/ 2 

-r-i t<-i 


and 


r 


3 

It | j e i(l-/W2 


a + i - r 


,1/2 J 


0<t< 1 


lim 77 (t + ie) = -s 
e-0+ 


i(l+/3)7r/2 


a + 1 - 1 


1/2 j 


-1 <t< 0 




It follows that 


lim arg tj (t + ie) = 
e-0 


(1 - /3)tt/ 2 0 < t < 1 

I(1+/3)tt/ 2 -1 < t < 0 



T?(t) 


i (1-/3)7 t/ 2+/3 tan" 1 V?-i| t> 1 


i (1-/3)7 t/2+j 3 tan" 1 ^-Vt 2 -T^ 


t < -1 


where the arc tangent must lie in the first and second quadrants. Also 


and 1 


(16a) 


(16b) 


(17a) 
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b(t) | = 1 t< -l, t> 1 


(17b) 


so that arg rj(t) goes from (1 - (3)ir/2 to jt/2 as t goes from 1 to °° along the real axis 
and from (1 + /3)^/2 to jt/2 as t goes from -1 to -°° along the real axis. Note, also, 
that when t = iX for X > 0, 



X > 0 


(18) 


Thus, 


arg T 7 (iX) =- X > 0 
2 


In fact, lim arg r] (t) = rr/2 everywhere in the upper-half plane, or 
t— OO 

lim ?j(t) = e i7r / 2 (19) 

t— «> 


Now consider £ . It is not known yet whether the boundary of the arc has a finite curva- 
ture at the stagnation point B or whether a sharp angle occurs at this point. There is 
no reason to exclude the latter case from consideration. Define 

lim <p(s^) - lim cp( s^) = /37T (20) 

sU on BC s^-~0 on AB 

that is, fin is the included angle between the two tangents to the arc surface at the stagna- 
tion point (see fig. 2). 

A classical result (ref. 4, p. 132) in conformal mapping is that in the vicinity of 
t = 0, dt/dz ~ iP ~ 1 . Because 


£ _ dW _ dW dt 
dz dt dz 

it follows, by using equation (12), that £ ~ iP in the vicinity of t = 0, which shows that, 
unless the arc has a cusp (/3 = 0) at the stagnation point, £ vanishes there as \P. It is 
also known from the nature of the problem that £ * 0, at any other point. Since the 
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velocity in the channel must become uniform far downstream of the arc, it follows from 
figures 2 and 4 that, in the upper-half t-plane, 

lim C(t) = V e i7r/2 (21) 

t-~°° ° 

Equation (15) shows that 77 (t) ~ tP in the vicinity of t = 0 and 77 (t) # 0, °° at any 
other point in the upper -half t-plane. Also, equation (19) holds as t — 

The interior of the first quadrant of the t-plane is denoted by T (i. e. , the first quad- 
rant of the t-plane excluding the positive real and imaginary axes), and the closure of T 
(that is, the first quadrant of the t-plane including the positive real an d im aginary axes) 
is dented by T. It is clear that half the flow (that is, the region in I BC HI of fig. 2) is 
mapped into T. 

From the preceding discussion, the function £ (t) /V Q ri (t) for t€ r is bounded away 
from 0 and <*> in T, is analytic in T, and approaches unity as t — °° for all t £ T. 
Introduce a new function of t, S2 = 9 + r (9 and r real) defined by 

?(t) = V 0 ??(t)e -in(t) (22) 

Then, S2 is analytic in T and bounded, and turns out to be continuous in T (ref. 5, 
Theorem 14. 19, p. 281 and Remarks 14. 20 C, p. 282). Also, since 

In 5ft) = -ifl(t) ter 

V Q r?(t) 

it follows from equations ( 20 ) and ( 21 ), that 

lim S2(t) = 0 teT (23) 

t— 


Also, figures 2 and 4 clearly show that £ must have the argument ir/2 on the positive 
imaginary axis of the t-plane. Since 

arg £ = arg 77 - 9 

it is clear from equation (19), 9 = 0 on the positive imaginary axis; that is, fi is a pure 
imaginary quantity on the positive imaginary axis. Then, by Schwarz's Principle of 
Reflection (ref. 4, pp. 44 or ref. 5, Theorem 11. 17, p. 230), S2(t) can be extended ana- 
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lytically to a function analytic in the entire upper-half plane (excluding the real axis) by 

fi(-t*) = -fi*(t) (24) 

Equations (15) and (24) show that £ (t) given by equation (22) then has the proper 
symmetry in the upper-half t-plane to describe the flow depicted in figure 2. Inasmuch 
as the free streamlines bounding the wake are mapped into the portion of the real axis 
outside the interval (-1, 1) (fig. 4), it can be concluded that [ £ (t) [ = V Q for t < -1 and 
t > 1. Equation (22) shows that 


|?(t)| = V 0 Mt)|e T(t) 


(25) 


thus, from equation (17a), 


e T ^ = 1 t< -1, t> 1 


which shows that 


r(t) = 0 t< -1, t> 1 (26) 

The complex conjugate velocity £ (t), given by equation (22) with the symmetry condition 
(eq. (24)) and the conditions of equations (23a) and (26), now satisfies all the specified 
velocity boundary conditions except those on the arc surface corresponding to the interval 
(-1, 1) in the t-plane. Equations (23) and (24) show that 

lim 0(t) = 0 (23) 

t— oo 


along all paths in the upper-half plane. 

Since S2(t) is analytic in the upper-half plane, vanishes at infinity, and is bounded on 
the real axis, fl(t) may be expressed everywhere in the upper-half plane in terms of its 
imaginary part on the real axis by using the Schwarz integral formula (ref. 3, p. 251), 
namely, 



JT 



ds 


14 



where 


s = Re (t) (27) 

By use of equation (26), the Schwarz formula becomes 

S?(t) = - f l ds (28) 

rt X\ s “ t 

It is easily seen that the function £2, given by equation (28) satisfies equations (23b) and 
(26) and that the condition 


t(-s) = t(s) (29) 

is necessary and sufficient for equation (24) to hold (which also guarantees that be a 
pure imaginary quantity on the positive imaginary axis). 

Observe that £(t), given by equations (22) and (28) with the symmetry condition 
(eq. (29)), satisfies all the boundary conditions on the velocity except the one on the 
boundary of the arc (-1, 1) in the t- plane. 

In order to find 0(s) for -1 < s ^ 1 from equation (28), the Plemelj formulas (also 
called Sohkosti formulas) for the limiting value of an analytic function (ref. 6, p. 25) are 
used. Corresponding to equation (28), the Plemelj formula appropriate to the upper-half 
t- plane is 

lim £2(t) = 0(s) + ir(s) 
t— s+ 


= ir(s) + P. V. — f 1 ) ds' -1 < s < 1 
^ J\ s' - s 

where t — s+ indicates that t approaches the real axis through values with positive 
imaginary parts, P. V. shows that the singular integral to the right is to be evaluated in 
the sense of a Cauchy Principal Value, and s' is a dummy variable of integration. Thus, 
0(s) is given on the interval (-1, 1) in terms of the values of r on this interval by 

0(s) = P. V. - f l ds’ -1 < s < 1 (30) 

* J x s' - s 
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Now, to complete the solution of the problem requires only the determination of r 
on the real interval (-1, 1), for which equation (6) is used. First, note that the element 
of length ds 1 is given by 


ds^ = |dz | 


where dz is to be taken along the real interval -1 < t < 1 in the t- plane, since s^ 
maps into this interval. Now, because 


dz^Mdt^r 1 ^ 

dW dt \dt / 


dt 


it follows that 


ds t = |? _1 (s)| 


dW 

dt 



-1 < s < 1 


In general, 


\r\s)\ = |?( S )| 


-i 


and, by definition, 


V = |C(s)| 


Using these results in equation (6) gives 


PV? dV o = jB cos cp 
s s 


dW 


dt 


ds 


t=s 


> -1 < s < 1 


v = |?(s)| 


(31) 


The surface of the arc is a streamline (see fig. 2); hence, 
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r 


<P 


0 < s< 1 


arg ? _1 (s) = 1 


:«-2 


■1 < s < 0 


but from equation (22), 


arg S"l(s) = 0 - arg 77(s) -1 < s < 1 


Using equation (16a) yields 


so that 


arg ? _1 (s) = 


0 - (1 - /3)ir/2 0 < s< 1 


0 - (1 + /3)7r/2 -1 < s < 0 


<P=< 


‘fl 


-0 


0 < s< 1 


■1 < s < 0 


(32) 


Integrating equation (31) from 0 to s (-1 < s < 1) and noting that V g is zero at s = 0 
(neglecting here the case of a cusp (j3 = 0) at the stagnation point) yield 


|?(s)| 3 h v 3 


3jB 

P 


/ 


cos cp( s’) 


dW 


dt 


Ids’ 


■1 < s < 1 


(33) 


t=s’ 


From equations (16) and (25), 
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|?(S)| = V 0 |t?(s) |e x ( s > 


P 


= V„ 


t(s) 


+ V7T> 


-1 < s < 1 


Substituting equations (12), (32), and (34) into equation (33) gives 


(34) 


3/3 


,3t(s) _ 


H-VlTs 2 / 


f COsLs.)^! - £ 
pv 2 X L 2 -l 1 + (s’) 2 /h 2 


ds' 0 < s < 1 


Setting s = 1 in equation (35) and using equation (26) show that 


-fzL. /•* «,(,*&)_ 

6jB ^o 4 V 2 'l + 


s 2 /h 2 


ds 


Eliminating 6jB^ 0 /pV^ between equations (35) and (36) leads to 


3/3 


»3t(s) _ 


a + ■ 


/ cos f^(s’) +— 1 — 

>_ * J l_+M 

/ cos/e +&) S_ 

1 2 ^ 1 * s 2 /l 


ds' 


ds 


From equations (29) and (30) or equivalently from equation (24), 


9(-s) = -9( s) 0 < s < 1 


(35) 


(36) 


0 < s < 1 (37) 


(38) 


so that 9 is an odd function of s on the interval (- 1, 1). In particular, 


0(0) = 0 


(39) 


18 


Equation (39) shows that in the neighborhood of s = 0 



2 

so that, in the neighborhood of s = 0, the right side of equation (37) behaves like s for 

3/3 

/3 * 1, whereas the left side behaves like s . It can be concluded from this that solu- 
tions to equation (37), which are continuous in the neighborhood of s = 0, exist only if 
/3 = 1 or /3 = 2/3. Thus, the arc surface is either a smooth curve (/3 = 1) in the neighbor- 
hood of the stagnation point, or else it is pointed with an included angle of 120° (compare 
with 120° wedge solution of ref. 2). Solving equation (37) for r yields 



0 < s < 1; /3 = 1, 2/3 (40) 

The problem is now essentially solved since the r given by equation (40) for positive s 
can be extended to negative values of s by the symmetry condition (eq. (29)). The re- 
sult can be substituted in equation (30). In this manner, a nonlinear integral equation for 
6 is obtained, and the solution to this equation can be obtained numerically. Equa- 
tions (40) and (29) can be used to calculate t(s) for -1 < s < -1. These values of r can 
then be substituted into equation (28) to determine 0 in the upper-half t-plane and, 
therefore, by equation (22) to determine the complex conjugate velocity £(t). Thus, 
knowing W from equation (13) and £ allows the calculation of all physical quantities, as 
explained in the discussion of equation (9) (see the Calculation of Stream Function sec- 
tion). The substitutions just described do not have to be carried out, since the numerical 
solutions to equations (29), (30), and (40) were obtained with the equations in the form in 
which they are now written. The solutions for an arc in an unbounded free stream can 


19 



be obtained by letting h — °° in equation (40) . 

Before the results are collected and summarized, the expressions for the most im- 
portant physical variables must be developed (the arc shape and the shape of the free 
streamlines), and the relations between the parameter h and the physical quantities of 
interest in the problem must be found. 


Boundaries of Arc and Wake 

Expressions for the shape of the bounding streamlines are obtained first. Notice that 

the surface of the arc BC in figure 2 is mapped into the interval 0 < s < 1 in the 

t- plane. Hence, the position of the arc surface, denoted by z = x + iy , is given by 

s s s 

equation (9) as 



0 < s < 1 


(41) 


where the origin of the coordinate system lies at the stagnation point. Substituting equa- 
tions (12), (16), and (22) into equation (41) yields 



and 
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This last expression can be simplified further. First, differentiate both sides of 
equation (35) with respect to s to get 



Using this expression to eliminate cos (0 + Pir/2) in equation (43) gives 



Denote the longest dimension of the arc in the direction parallel to the flow by d; thus, 

y = -d at s = 1. From equation (26), t( 1) = 0, so that equation (44) shows 
s 



2jB 


or 
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( 45 ) 


which is the equilibrium condition satisfied by the parameters of the magnetically bal- 
anced arc. Using equations (45) and (36) in equation (42) finally leads to 


/ 

fs = 2 ^ 

d 3 


cos [0 - (1 - P)n/2]e~ 


/■ •-(•*?): 


v i + Yi - s* 
1 + s 2 /h 2 


+ s 2 /h 2 




s 1_/3 ds 


0 < s < 1 (46a) 


Z_ = 2 
d 3 


cos [0 - (1 - P)n/2]e~ 


/ •“(•*?): 


^1 + T1 - s- 
1 + s 2 /h 2 


+ s 2 /h 2 


s 1-/3 ds 


(46b) 


and rewriting equation (44) in terms of d yields 


.1 + T 1 - s' 


a 2r(s) 


0 < s — 1 


Having obtained a set of parametric equations for the bounding surface of the arc, the 
next step is to obtain an expression for the boundary of the wake. Notice that the bound- 
ing streamline of the wake CH in figure 2 is mapped into the interval s> 1 in the 
t-plane. Hence, the position of the wake surface defined as = x^ + iy^ is given by 
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z f - (z s> 


+ r rW— ) 

s=i X uy t= s 


ds s > 1 


(48) 


Substituting equations (12), (17), (22), and (26) into equation (48) leads to 

-/3)7r/2-/3 tan" 


5 f ' 



i 0-(l-/3)7r/2-/3 tan -1 Vs 2 -l 


i 2 /u 2 

1 + s /h 


ds s 2: 1 


The real and imaginary parts of are 


x f = l + 


^ f cos [ 

v o X L 


0 - (1 - /3)7r/2 - /3 tan' 




1 + S 2 /h 2 


and 


V (49) 


2 ^o /* s r 

y f = -d + J sin - (1 


- f3)ir/2 - /3 tan' 


_ 1 + 


2 /, 2 
s /h 


ds 


where Z is given by equation (46b). Using equations (45) and (36) in equation (49) gives 




6 - (1 - /3)n/2 - 2 tan" 




d d 3 


ds 


J 1 + s 2 /h 2 


/“<••? >77 


s > 1 (50) 


2 a2 

s /h 


ds 
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Ill" 


l!~ i + * 

d 3 


/ 


sin 0 - (1 - P)n/2 - )3 tan’ 1 f s 2 - 1 


cos e +&) 5- 


1 + s 2 /h 2 


/ 


s > 1 (51) 


2 'l + s 2 /h 2 


Equations (46a) and (47) and equations (50) and (51) give the shape of the arc and wake, 
respectively. Taking the limit h — °° in these equations gives the results for an arc in 
an unbounded free stream. 


Physical Parameters 

To obtain an expression for V 00 /V 0 in terms of h, note that far upstream of the arc 
£ = iV . In the t- plane, this point corresponds to t = ih, and so equation (22) shows that 


iV^ = V 0 r?(ih)e 


ifi(ih) 


but from equation (28), 


S2(ih) = - 


1 r 

• j i 


s - ih 


!/ 


_Sll2L ds + i!! 

s 2 + h 2 ' 


/ 


2 h 2 
s + h 


2 2 

Now equation (29) shows that r(s)/(s + h ) is an even function of s, and thus 
2 2 

st(s)/(s + h ) is an odd function of s. The first integral on the right of equation (53) is, 
therefore, zero. Introducing equations (18) and (53) into equation (52) yields 


h + 1 + 




l+s 2 /h 2 


■ 


ii i ii m 



Equation (54) shows that in the limit h — °°, V Q = V^. Thus, in the case of an unbounded 
arc, the velocity far upstream must be equal to the velocity on the boundary streamline 
of the wake (or the velocity far downstream), which could have been anticipated before- 
hand from known results about unbounded flows. The characteristic dimension of the arc 
d can be related to the channel width D by first eliminating 4 > q between equations (14) 
and (36) to yield 


f\os(e^)- 

3jBv - D X ' 2 y i+ 


s 2 /h 2 


ds 


Using equation (45) gives 


d Zo 

D V 


2nh 




s 2 /h 2 


ds 


or, with V 00 /V 0 eliminated by equation (54), 



which again shows that D — <» when h — 

The cross-sectional area A of the arc is given by 



Substituting equations (46a) and (47) into the expression for A results in 
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P s'-P 


A=id 2 


cos [0 - (1 - /3)7 t/ 2] e T ^1 - Vl - S 2 ^ — 


s 2 /h 2 


ds 


/ cos ( e+ f)- 


s 2 /h 2 


ds 


and equation (46b) is used to get 


0 o l-0 


A 
22 d 


/ cos [0 - (1 - 0)ff/2]e T /l - Vl - s 2 j — - 

V 1 + 


2 ,.2 
s /h 


ds 


/ 


cos [0 - (1 - 0)7r/2]e' 


7 


0 1-0 


+ s 2 /h 2 


ds 


Combining equation (56) with equations (45) and (54) yields the drag coefficient 


SB 


(H) w 


VF77 + iV . Th S .-W 


r t(s) 

h l+s 2 /h 2 


ds 


0 s i-0 


/ cos [0 - (1 - 0)7T/2]e T ^1 - Vl - s 2 j 

1 + 


s 2 /h 2 


0 , 1-0 


cos [0 - (1 - 0)7r/2]e" T ^1 + Vl - s 2 ^ — - 


s 2 /h 2 


ds 


ds 


(56) 


(57) 


where S= jA is the total current in the arc. 


Summary of Calculated Quantities 


For convenience, the results are now collected. If 0 = 1, the boundary of the arc is 
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smooth at the stagnation point; if /3 = 2/3, the boundary of the arc is pointed with 120° 
included angle between the sides of the arc at the stagnation point. The parameter /3 
cannot have any other values. The following is a summary of the final results: 


t(s) = /3 In 


LVTTTli 


In 


, /*'“("*?) 77 


s 2 /h 2 


ds 


,3/3 


/“■(•♦f)— Vi 

4 V 2 ' 1 + s 2 /h 2 


ds 



0 < s < 1; (3=1, 2/3 

(40) 

il 


(29) 

/ T < s ’> ds’ 
' \ s’ - s 

-1 < s < 1 

(30) 


Equations (45) and (54) yield 
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cos [0 - (1 - /3)7r/2]e 




l_ = 2 
d 3 


/ cos ( e+ t)- 


1 + s 2 /h 2 


ds 


s 2 /h 2 


(46b) 


Taking the real part of equation (28) results in 


0(s)=- 


= 1 r 1 r(s') 

7T /i S' - S 


ds’ s > 1, s < -1 


(59) 


_s _ 2 
d 3 



cos [0 - (1 - P)it/2]e 




/“(’•St: 


i 2 /, 2 

1 + s /h 


ds 


s x K ds 


0 < s < 1 (46a) 


s 2 /h 2 


t 2/S 


,2t(s) 


0 < s < 1 


d \l + VT7s 2 > 


(47) 


= Z_ + 2 
d d 3 


/“[■ 


0 - (1 - /3)tt/ 2 - /3 tan’ 




J 1 + s 2 /h 2 


ds 


/ 


s — 1 (50) 


cos 1 0 + ^ | ds 


2 / 1 + s 2 /h 2 
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Finally, the drag coefficient yB/(pV^ 0 /2)(2l) is 



NUMERICAL PROCEDURE 

Equations ( 30 ) and ( 40 ) were solved numerically with |3 = 2/3 for various values of 
h (including h = °°) and with /3 = 1 for h = <*>. These equations were solved by iteration. 
First, 0 was guessed as 0(s) « s and substituted into equation ( 40 ) to calculate r(s). 
This t(s) was used in equation ( 30 ) to calculate a new 0(s) which was then used in equa- 
tion ( 40 ) to calculate new values of r(s). The procedure was repeated until convergence 
was obtained. In the application of this iteration procedure, certain special techniques 
had to be used. First, the numerical integration scheme used with equation ( 30 ) had to 
be tailored to take the principal value of the integral; that is, it had to be set up so that 
the singularities on both sides of the point s f = s cancelled one another in the principal- 
value sense. Second, if any small bumps developed in the curve t against s during the 
course of the iteration process, it was necessary to smooth these out before the values of 
r were substituted into equation ( 30 ) to calculate new values of 0. This smoothing out 
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was done because the singularity in equation (30) would tend to make these small bumps 
grow and thus lead to an unstable iteration process. Finally, it was necessary to use an 
averaged iteration procedure with respect to a weight factor e (ref. 4, p. 216); that is, 
instead of equation (30) being used directly in the iteration processes, it was necessary 
to rewrite it in the form 


0(s) = e0(s) + (1 - e)P. V. - Z* 1 ds' 

* JO s ’ - s 

so that, if 6^ and are the values of 6 and r obtained in the n** 1 step of the iter- 
ation process, the n + 1 step value of 6 is obtained from 

g ( n +!) ( s ) = e0 (n)( s ) + (l _ e )p. v. A C T ^( s> ) 

* J-\ S S 

When these techniques were used, the iteration process converged to three significant 
figures in less than 10 iterations with e = 0. 5. 

RESULTS AND DISCUSSION 

The analysis shows that there are two possible solutions to the problem of determin- 
ing the shape of an arc with the model used for the internal structure. One of these arcs 
has a sharp leading edge (figs. 6 and 7) with a 120° included angle, and the other has a 
blunt leading edge (fig. 8). The solution for the arc with the sharp leading edge (with the 
channel walls separated by six or more arc widths, which corresponds to the experimen- 
tal conditions) predicts an arc width-thickness (in the direction of flow) ratio of about 2 
and a value of the experimentally observed parameter JB/(p\^ 0 /2)(2l) between 0. 6 and 
0. 8 (fig. 9). These results are in agreement with the experimental results of Roman 
and Myers (except those results at the lower velocities) (ref. 1) and the results of Myers 
(ref. 7). The dimensionless parameters shown in figure 9 are uniquely related to the dis- 
tance h in the t-plane (fig. 4) according to the relations listed in the section Summary 
of Calculated Quantities. In particular, the parameter h determines the shape of the 
arc. Even for an infinite stream (h = °°), a change in the stream velocity alone, for ex- 
ample, does not alter the shape of the arc because the magnetic field or current must be 
adjusted to maintain stability. 

The arc with the blunt leading edge (/3 = 1) had too large a width-thickness ratio to 
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Figure 6. - Shape of arc with sharp leading edge in infinite stream. Distance between 
stagnation point and point at upstream infinity in t-plane, <»; channel width, 
included angle between two tangents to arc surface at stagnation point divided by 7r, 
2/3. 



y/d 


Figure 7. - Shape of arc and wake with sharp leading edge. Distance between stagnation 
point and point at upstream infinity in t-plane, 1; channel width, 5.88; included 
angle between two tangents to arc surface at stagnation point divided by n, 2/3. 
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Figure 8. - Shape of arc with blunt leading edge in infinite stream. Distance between stagnation point and point at 
upstream infinity in t-plane, «>; channel width, included angle between two tangents to arc surface at 
stagnation point divided by 7r, 1. 
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Figure 9. - Width -thickness ratio and drag coefficient for arc with sharp leading 
edge as function of ratio of channel width to arc width. Included angle between 
two tangents to arc surface at stagnation point divided by re, 2/3. 
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Figure 10. - Width-thickness ratio and dimensionless 
ponderomotive force as function of channel size for 
arc with blunt leading edge. Included angle between 
two tangents to arc surface at stagnation point 
divided by 7r, 1. 


agree with the experiments. The geometric and dynamic parameters for this type of arc 
are shown in figure 10. The actual bluntness observed in experiments must be attributed 
to transport effects such as viscosity. The model also predicts the broad wake region 
observed experimentally. At this stage, it is not as yet possible to make any more com- 
parisons with the experimental results. 

The techniques developed herein can be used easily with more general models of the 
internal structure of the arc to obtain more accurate results. 


CONCLUSION 

A previous analysis of the shape of a magnetically balanced arc was improved by 
giving a more accurate description of the external flow field. Two possible solutions 
were obtained for the model used. One solution was in good agreement with experimental 
results whereas the other was not. The model used for the internal structure of the arc 
is expected to be most accurate at the higher flow velocities. 
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The techniques developed in this report can be used easily with more general models 
of the internal structure of the arc to obtain more accurate results. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 31, 1968, 

129-01-11-05-22. 
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